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RESUMEN 



EI12 de agosto de 2014, se reglstro un sismo de magnltud 5.1 , a una profundldad focal de 4 km., 
en el segmento de falla Bellavlsta Catequllla, el mismo que fue reglstrado en 8 estaclones locallzadas en 
la cludad de Quito. Estas se encuentran ubicadas a distanclas epicentrales entre 12 y 1 9 km. 

En este artfculo se comparan las aceleraclones maximas obtenldas en campo llbre, con las que 
se obtlenen al emplear las ecuaclones de movlmlentos fuertes de Campbell y Borzognia (2013) y el de 
Zhao et al. (2006). Para ello prevlamente se determlna un piano de ruptura del sIsmo, utillzando las 
ecuaclones propuestas por Leonard (2010) y la geometria de las fallas degas propuestas por Alvarado 
(2014) 



ABSTRACT 



On August 12"^ 2014, a magnitude 5.1 earthquake occurred at a depth of 4 km. In the Bellavlsta 
Catequllla fault segment. This event was recorded by 8 strong-motion stations located between 12 and 19 
km from the epicenter. In the city of Quito. 

In this article, the maximum ground accelerations recorded In free field are compared with the 
accelerations estimated using the models by Campbell y Borzognia (2013) and Zhao etal. (2006). To this 
end, the earthquake rupture plane Is determined using the equations proposed by Leonard (2010) and the 
geometry of the blind fault system of Quito proposed by Alvarado (2014). 
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1. INTRODUCCION 



El martes 12 de agosto de 2014, la naturaleza volvio a recordar a los habitantes de la ciudad de 
Quito, que sus viviendas se encuentran sobre las fallas ciegas de la ciudad o muy proximas a ellas, al 
ocurrir un sismo de magnitud 5.1, cuyo epicentro, preliminarmente fue ubicado en las siguientes 
coordenadas: Latitud -0.06 S; Longitud -78.38 W., y a una profundidad focal de 4 km., de acuerdo al 
Institute Geofisico de la Politecnica Nacional, IPGN. En una revision posterior el epicentro se ubico con 
Latitud 0.05 S; Longitud -78.42 W y una profundidad focal de 5 km. 

En la figura 1 , se presentan varias ubicaciones del epicentro proporcionadas por el ISC 
(International Seismological Center, Agenda Sismologica GFZ de Alemania); USGS (Servicio Geologico 
de los Estados Unidos); NEIC (National Earthquake Information Center) y MOS (Agenda sismologica de 
Rusia). Los mayores danos se registraron en la zona propuesta por el IGPN, para el epicentro del sismo. 

Los sismos asociados a estas fallas inversas son de tipo impulsive y tienen las siguientes 
caracteristicas: son de corta duracion; causan demasiado dano en la zona epicentral cuando la magnitud 
es mayor que 6, pero se atenuan muy rapidamente con la distancia. Para tener una idea del poder 
destructive de los sismos impulsivos basta mencionar los sismos de Northridge de 1 994 (Estados Unidos) 
y el de Kobe de 1 995 (Japon). 



^JU'J III 




78°40'W 78°35'W 78°30'W 78°25'W 78°20'W 78°15'W 78°10'W 78°5'W 

Figura 1 Localizacion del epicentro del sismo del 12 de agosto de 2014, segun distintas 

agendas. 

Con esta acotacion, en la tabia 1, se indican las aceraciones maximas registradas del sismo del 
12 de agosto de 2014, a mas de 10 km., de distancia epicentral. La mayor aceleracion es de 48 gal y fue 
registrada en la estacion PRAM, que se encuentra muy proxima al antiguo aeropuerto de Quito. No se 
disponen de registros en la zona epicentral pero los correspondientes valores de PGA debieron ser mas 
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grandes que las indicadas en la tabia 1 , que son bastante bajas y como era de esperarse solo causaron 
un gran susto en la mayor parte de la poblacion de Quito. 



TabIa 1 Valores de PGA y ubicacion de estaciones (Fuente: IGEPN) 



Estacion 


Lat_estac 


Long_estac 


PGA 

(gal) 


Dist 
(Km) 


PRAM 


-0,14 


-78,49 


48 


13 


FENY 


-0,15 


-78,48 


38 


12 


24MA 


-0,18 


-78,47 


18 


15 


lESS 


-0,18 


-78,49 


31 


16 


IRDE 


-0,20 


-78,48 


20 


17 


CIRC 


-0,20 


-78,49 


41 


17 


GMEA 


-0,20 


-78,49 


31 


17 


AEPN 


-0,21 


'78,49 


21 


19 



En los valles de los Chillos y de Tumbaco, que se encuentran aproximadamente a 30 km. del 
epicentre se sintio como un pequeno temblor. Esto se refleja en el mapa de isosistas presentado en la 
figura 2 donde puede apreciarse en la zona epicentral su intensidad fue 6 en la escala macro sismica 
europea, EMS-98 mientras que en la mayor parte de la ciudad fue de 4, oscilando en los valles entre 3 y 
4. 




-79.5 -79 -78.5 -78 -77.5 

Intensidad (EMS-98) 




3 4 5 



Figura 2 Mapa de isosistas del sismo del 12 de agosto de 2014. Fuente: Singaucho, 2014. 

El sismo del 12 de agosto de 2014, puede ser considerado como un recordatorio de que la ciudad 
de Quito se asienta sobre un sistema de fallas inversas que en la mayor parte no tiene afloramiento 
superficial, pero hay una serie de evidencias que demuestran su existencia, como es la forma en que se 
encuentra el material volcanico a la salida del tunel de Guayasamin, en direccion al valle de Tumbaco, o 
la alineacion de las colinas de llumbisi y Puengasi. 
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2. REPLICAS DEL SISMO Y NECESIDAD DEL ESTUDIO 



Mas de 60 replicas del sismo de Quito de Agosto de 2014 han sido registradas. Las de mayor 
magnitud se indican en la tabia 2. La ubicacion de los epicentros del evento principal y las replicas se 
muestra en la figura 3. La cantidad indicada corresponde a la profundidad focal. En la figura 3, se aprecian 
cuatro de los cinco segmentos que conforman las fallas ciegas de Quito y de sur a norte son: Puengasi; 
llumbisi-La Bota (ILB); Carcelen-EI Inca (CEI); Bellavista-Catequilla (BC). Un poco mas al norte se halla el 
segmento de Tangahuilla. 

Ta bIa 2 Ubicacion de las replicas de mayor magnitud del sismo del 12 de agos to de 2014 



Fecna 


Coordenadas 


iviagniiuu 


Profundidad 
Focal 


Latitud 


Longitud 


12/08/2014 


-0,05 


-78,42 


5,1 


5,0 


12/08/2014 


-0,05 


-78,4 


4,0 


5,0 


12/08/2014 


-0,06 


-78,42 


3,7 


5,0 


12/08/2014 


-0,05 


-78,43 


3,3 


6,0 


12/08/2014 


-0,05 


-78,41 


3,1 


6,0 


13/08/2014 


-0,05 


-78,41 


3,4 


4,0 


13/08/2014 


-0,05 


-78,42 


3,3 


2,0 


13/08/2014 


-0,06 


-78,42 


3,9 


5,0 


13/08/2014 


-0,05 


-78,42 


3,8 


4,0 


15/08/2014 


-0,07 


-78,44 


3,0 


11,0 


16/08/2014 


-0,05 


-78,41 


4,8 


4,0 


16/08/2014 


-0,06 


-78,41 


3,1 


6,0 


16/08/2014 


-0,05 


-78,42 


3,7 


7,0 


16/08/2014 


-0,05 


-78,41 


3,1 


5,0 


16/08/2014 


-0,04 


-78,41 


3,1 


5,0 


17/08/2014 


-0,03 


-78,43 


3,0 


5,0 


17/08/2014 


-0,07 


-78,43 


3,3 


7,0 


17/08/2014 


-0,05 


-78,42 


4,0 


6,0 


17/08/2014 


-0,07 


-78,42 


3,0 


5,0 


21/08/2014 


-0,07 


-78,41 


3,0 


4,0 


22/08/2014 


-0,08 


-78,43 


3,1 


7,0 


25/08/2014 


-0,08 


-78,43 


3,3 


10,0 



En base a esta informacion se obtuvo la ecuacion de recurrencia propuesta por Gutenberg y 
Richter, la misma que se presenta en la figura 4, donde se aprecia que existe una muy buena correlacion 
entre la magnitud y el numero de sismos. La ecuacion de ajuste encontrada es la siguiente. 

log(Af) = 3.0355 - 0.5935 M 

Donde N, es el numero de sismos de magnitud igual o mayor que M. 

En la figura 3 se aprecia que este sismo esta asociado al segmento: Bellavista-Catequilla. Mas 
adelante se vera que el 10 de agosto de 1990 se produjo un sismo de magnitud M„ = 5.3 a pocos 
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kilometros del epicentre del sismo de 2014, en el segmento de falla BC. Por lo tanto, en 24 anos se ha 
vuelto a liberar energia. 



0*05'0"N 



O'O'O" 




5.of n'e.o 

11.0» • • •5.0 CATEQUILLA-BELLAVISTA 
5:0 4.0 

r 




Figura 3 Localizacion del evento principal y replicas del sismo del 12 de agosto de 2014 



Alvarado ef a/. (2014) en base al modelo del Terremoto Caracterfstico, determina para el segmento 
BC un terremoto de magnitud 6.3, con un periodo de recurrencia de 1 83 anos. Este calculo es inquietante 
para los Proyectistas Estructurales, ya que el periodo de recurrencia es pequeno y la magnitud es alta para 
un probable sismo superficial que puede producirse dado que el espectro que genera un sismo de estas 
caracteristicas es alto y el periodo de recurrencia pequeno, se debe verificar que la estructura tenga muy 
poca incursion en el range no lineal. Aguiar (2003). 
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En este articulo se presenta el modelo de Terremoto Caracterfstico para determinar el perfodo de 
recurrencia para el segmento de falla BC y demostrar que efectivamente el periodo de recurrencia con 
este modelo es de 183 anos, pero es muy importante conocer las hipotesis del modelo para no alarmarse 
con las magnitudes esperadas en los segmentos de las fallas ciegas de Quito, asociadas a periodos de 
recurrencia cortos. 

Sin embargo es necesario verificar el desempeno estructural ante un sismo muy especifico, en la 
ciudad de Quito, ya que el 70% de sus construcciones se encuentran sobre o muy cerca de las fallas 
ciegas. Claro esta que el diseno de la estructura debe realizarse para los espectros encontrados en el 
estudio de la microzonificacion sismica de Quito. ERN (2012), Aguiar (2013). 



1,6 n 




2,5 3,5 4,5 5,5 



Figura 4 Ajuste de las replicas con el modelo de Gutenberg y Richter. 
3. TECTONICA LOCAL 

Alvarado et al. (2014), al trabajar con las ecuaciones de Leonard (2010) determinan dos 
magnitudes para cada uno de los segmentos de las fallas ciegas inversas de Quito, las mismas que se 
indican en la tabia 3. En base al modelo del Terremoto Caracterfstico (Wesnousky, 1986) encuentran el 
perfodo de recurrencia. 



TabIa 3 Segmentos de las fallas ciegas de Quito. Alvarado et al. (2014). 



Segmento 


Area 


Magnitud 


Longitud de 


Magnitud 


Perfodo de 




ruptura 


(RA) 


Superficie 


(SRL) 


Recurrencia 




(km)2 




de Ruptura 




(anos) 








(km) 






Puengasi 


259 


6.4 


22 


6.4 


188 


ILB 


176 


6.2 


15 


6.2 


138 


CEI 


82 


5.9 


7 


5.7 


105 


BC 


191 


6.3 


17.5 


6.3 


183 


Tangahuilla 


108 


6.0 


12 


6.0 


115 



Las ecuaciones propuestas por Leonard (2010) para fiallar la magnitud en funcion de la longitud 
de ruptura L, con la que se encontro la magnitud (SRL) y para encontrar la magnitud en base al area de 
ruptura A, con la que determine la magnitud (RA) son las siguientes, para el caso de levantamiento de un 
bloque. 



= 1.52 Iog(L) + 4.4 



(1) 
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M„ = 1.01ogU) + 4 



(2) 



En el presente caso, el bloque de Quito se esta levantando con respecto a los Valles de los Chillos 
y Tumbaco por compresion tectonica. 
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Figura 5 Zonas fuentes locales y replicas del sismo del 12 de agosto de 2014. 



En base al area y longitud de ruptura, indicados en la tabia 3 y en base al angulo de buzamiento 
de 55°, Alvarado et al. (2014) encontro zonas fuentes locales, para cada uno de los segmentos, de las 
fallas ciegas de Quito. En la figura 5 se presenta la proyeccion en la superficie de los pianos de falla. De 
sur a norte los pianos corresponden a los segmentos de Puengasi, llumbisi la Bota, Carcelen el Inca y 
Bellavista Catequilla. En esta ultima zona fuente local se aprecia que el epicentro del sismo del 12 de 
agosto de 2014 y la mayor parte de las replicas se ubican dentro de la zona fuente BC. 



4. 



MODELO DEL TERREMOTO CARACTERISTICO 



El modelo del Terremoto Caracteristico, TC supone que en la falla se producen rupturas que 
generan un sismo de la misma magnitud con una periodicidad perfecta como se muestra en la figura 6. Se 
puede pensar que se rompe toda la falla o un segmento pero la cantidad que se considere se va a volver 
a reproducir en un tiempo constante (Wesnusky 1 986). 
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A modo de ejempio en la figura 6 se ha supuesto que cada 50 anos se genera un sismo que tiene 
una magnitud maxima Mjc (Magnitud Caracteristica). Luego de un terremoto de magnitud Mjc la escala 
de tiempo vuelve a cero y empieza a acumularse la energfa. Dado que la falla tiene una tasa de movimiento 
constante (velocidad) u, el nuevo sismo de magnitud Mjc se producira en 50 anos, para el caso del ejempio 
de la figura 6. Es decir en el ejempio mostrado el periodo de recurrencia es de 50 anos. Si en este lapso 
de tiempo se registra un sismo de menor magnitud, esta situacion no modifica el periodo de retorno con el 
modelo de Terremoto Caracterfstico. 




0 50 100 150 200 

Tiempo (anos) 



Figura 6 Modelo del Terremoto Caracterfstico. 

En el caso del segmento de falla BC este sismo tiene una magnitud de 6.3, considerando toda el 
area de ruptura, y el periodo de recurrencia asociado a esta magnitud es de 183 arios. El modelo supone 
que una vez que se da el sismo de 6.3, se debe esperar 183 arios para tener otro sismo de la misma 
magnitud; evidentemente si se han registrado ya dos sismos en BC, ocurridos en 1990 y 2014 es logico 
pensar que el periodo de recurrencia va a ser mayor de los 183 anos, porque ha habido una importante 
liberacion de energfa. 

El modelo de Terremoto Caracterfstico iguala el momento sfsmico acumulado en una falla 
Moacumuiado durante un tiempo de recurrencia medio T con el momento sfsmico liberado Moiiberado en el 
terremoto de magnitud caracterfstica. Lo que significa que toda la energfa acumulada en el tiempo de 
recurrencia T se libera en el sismo de magnitud caracterfstica. Numericamente se tiene. 

Mo(.MOacumulado) = "(Mrc) * Mq {MOaberado) (3) 

Donde Mq{Mo acumulado) es el momento total acumulado en la falla; n(Mj,c) es la tasa de momento 
sfsmico anual; Mq (Moiii,ej.ado) es el momento sfsmico liberado. Por otra parte el perfodo de recurrencia no 
es mas que la inversa de la tasa de momento sfsmico anual. 

Ahora bien, a partir de la longitud o del area de ruptura, se puede encontrar la magnitud maxima 
esperada en una fuente. En la tabia 3 se ha encontrado esta magnitud para cada uno de los segmentos 
de las fallas ciegas de Quito empleando las ecuaciones de Leonard (2010). Se pudo trabajar con otras 
relaciones como las propuestas por Wells y Coppersmith (1994); String et al. 2002, entre otras. 

Por medio de la expresion de Hanks y Kanamori (1979), se puede conocer el momento sfsmico 
(co-sfsmico) que se liberarfa en el terremoto de momento maximo a partir del dato de magnitud My^. 

= ^log(Mo) - 10.7 (4) 

Con los datos de la tabIa 3, para el segmento BC se tiene = 6.3, con lo que el momento sfsmico 
liberado, es. 

(Mw+10.7)3 

Mo = 10 2 = 3.1623 * lO^s dyn cm = 3.1623 * W^^ N m 
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Por otro lado, con la hipotesis conservadora de que en todo el piano de falla se esta acumulando 
la energia en forma uniforme, el momento total acumulado de acuerdo al modelo de Brune (1968) es. 

Mo = fiuA (5) 

Donde n es el modulo de deformacion de las rocas durante el sismo: u es el desplazam lento de la 
falla; A es el area de ruptura. En la tabia 3 se considero que (i = 3* 10^° Pa que es igual a ju = 3 * lo^ 
Al derivar la ecuacion (5) con respecto al tiempo se encuentra. 

Mo = yLuA jgj 



Donde Mo es la tasa de momento sismico anual; u es la tasa de movimiento de la falla, que para 
del segmento BC se considero li = 3 — = 0.003 — . 

ano ano 

El area de ruptura A = 191 fem^ = 191 * 10^ m^. Al reemplazar los datos en la ecuacion (6) se 



halla. 



Mo = 3* 10^ * 0.003 * 191 * 10^ = 1.719 * 10^^ — = 1.6846 * 10^^ 

ano ano 

Finalmente el periodo de recurrencia T se halla con la siguiente ecuacion: 

Mo 

T = ^ (7) 
Mo ^ ' 

Donde Mo es el momento sismico liberado durante el sismo y Mo es la tasa de acumulacion del 
momento sismico. 

3.1623 * 10^^ Nm 
T = Tj — = 187.7 anos 

1.6846 * 1016 

ano 

El valor estimado es del orden del calculado por Alvarado et al. (201 4) indicado en la tabIa 3. Este 
es de 183 anos, siendo la diferencia con el estudiado de menos de 5 anos, atribuible a los decimales de 
calculo. 



5. AREA DEL PLANO DE RUPTURA 



Para la determinacion del area del piano de ruptura del sismo del 12 de agosto de 2014, se utiiiza 
las ecuaciones propuestas por Leonard (2010) para el caso de un sismo asociado a una componente 
vertical del movimiento de suelo; estas ecuaciones son la (1) y la (2) indicadas al principio. 

Al remplazar M^ = 5.1 en las ecuaciones 1 y 2 se halla que la longitud de ruptura es de 2.9 km, y 

el area de ruptura es 12.6 km^. Al modelar el area del piano de ruptura como un rectangulo se encuentra 
que el ancho del piano de ruptura es W = 4.4 km. En la figura 7 se presenta la proyeccion del piano de 
ruptura en la superficie y la ubicacion del sismo principal y las replicas, segun la ubicacion dada por el 
IGPN. 
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-0°0" 



Figura 7 Piano de ruptura del sismo del 12 de agosto de 2014 



6. MODELOS DE MOVIMIENTO FUERTE 
6.1 . Modelo de Campbell y Borzognia 2013 

En 2008, Campbell y Borzognia, propusieron la siguiente ecuacion para determinar la aceleracion 
del suelo: 

In = fynag + fdis + fflt + fhng + fsite + fsed (8) 

Donde f^ag es una funcion de la magnitud M^^; f^i^ es una funcion de la magnitud y de la distancia 
del sitio al piano de ruptura Rmp', ffu depende de la magnitud y del tipo de falla; f^ng es una funcion que 
depende del efecto hanginwall, la ciudad de Quito se encuentra en el bloque que se esta levantando; fsue 
es una funcion que toma en cuenta los efectos de sitio y depende fundamentalmente de V^^o velocidad de 
la onda de corte en los primeros 30 metros; fsed es una funcion que toma en cuenta el efecto de cuenca, 
depende de Z2 5 que es la velocidad de la onda de corte a una profundidad de 2.5 km; Yes la aceleracion 
del suelo esperada como una fraccion de la gravedad. 

Ahora, bien, la nueva ecuacion propuesta por C & B 2013 es la siguiente. 

In V = /maa + fdis + fflt + fhng + fsite + fsed + fhyp + fdip + fdtn (9) 

El significado de los terminos que faltan es como sigue: f^yp es una funcion que depende de la 
profundidad focal Z^y^; f^ij, es una funcion que depende del angulo de buzamiento de la falla 5 y de la 
magnitud; f^tn es una funcion de la atenuacion anelastica. 
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Como se podra apreciar la ecuacion propuesta en 2013, es un complemento de la ecuacion 
propuesta en 2008 y ademas se ha incorporado a la base de dates los siguientes sismos ocurridos en: 
2006, California M=6; 2008, China M=7.9; 2009, Italia M=6.3; 2010, Mexico M=7.0; 2010, Nueva Zelanda 
M=7.0; 2011 Nueva Zelanda M=6.2. Logicamente toda esta informacion mas el nuevo conocimiento 
cientifico que se da dia a dia, llevo a C & B a presentar una nueva ecuacion con la cual se pueden 
determinar espectros de respuesta elastica asociados a un determinada falla geologica cortical o leyes de 
atenuacion de movimiento del suelo para diferentes periodos. 



6.2. Modelo de Zhao ef al. 2006 



El modelo de Zhao et al. 2006, permite predecir el movimiento del suelo y encontrar espectros para 
tres regimes tectonicos: corteza activa, subduccion interplaca y suduccion intraplaca. La mayor parte de la 
base de dates de este modelo proviene de estas dos ultimas fuentes sismicas y en su mayor parte fueron 
registrados en Japon. La ecuacion de este modelo es la siguiente: 

ln(Y) = aM^ + bX- ln(r) + e (ft - hc)S„ + + S, + Ss + Ssl In^) + Q (10) 

Donde es la magnitud; X es la distancia a la fuente en km., que en el anterior modelo se 
denomina Ry^p', r es una funcion que depende de la magnitud y de la distancia; h es la profundidad focal 
en km.; es un factor que solo se considera en sismos corticales provenientes de falla inversa; S, factor 
de correccion per distancia para sismos de interface; Ss,Ssl factores de correccion per distancia para 
sismos intraplaca; Q es un factor que considera los efectos de sitio per efecto del suelo; la clasificacion de 
suelos que considera el modelo es la de la guia NEHRP en que se clasifica el suelo en funcion de la 
velocidad de la onda de corte Vs^q. Los restantes parametros de la ecuacion (10) vienen tabulados en 
funcion del periodo. 

Tanto con las ecuaciones indicadas, para los modelos de C & B (2013) como la de Zhao et al. 
(2006) no se ha incluido las funciones con las cuales se evaluan la dispersion que tiene la base de dates 
utilizadas por lo que la probabilidad de excedencia al usar las ecuaciones en la forma presentada es del 
50%. 



7. CALCULO DE ACELERACIONES MAXIMAS 



En la tabia 4, se presentan las variables que se requieren para calcular la aceleracion maxima del 
suelo con los modelos indicados. Falta por indicar el significado de es la distancia mas cercana a la 
proyeccion de la superficie del borde superior del piano de ruptura medido perpendicularmente a su rumbo 
promedio; Ry^ es la distancia minima a la proyeccion de la superficie del piano de ruptura (es la distancia 
Joyner-Boore). 

El calculo se hizo considerando una magnitud M^v = 5.1; y profundidad focal Z^y^ = 5 km. La 
profundidad del techo del piano de ruptura co-sismica Z^or = 3.2 km., la profundidad al fondo del piano de 
ruptura Z^ot = 6.78 km. 

Lamentablemente no se conoce el valor de Vs3o para los lugares en que se hallan las estaciones, 
pero en la microzonificacion sismica de Quito, se tienen valores para sitios cercanos por lo que se decidio 
trabajar con un solo valor de ¥^30 = 400 m/s, en el modelo de C & B y corresponde a un suelo tipo C, en 
el modelo de Zhao et al. 2006. En la tabIa 5 se presentan las aceleraciones registradas y las calculadas 
con los dos modelos. 
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Tabia 4 Valores de Rmv.Rx'^jb' Pue se necesitan en los modelos 



Estacion 










(km ) 


(km.) 


(km ) 


PRAM 


14.0135 


13.1786 


1 1 .4452 


FENY 


13.2830 


12.3246 


10.7952 


24MA 


15.3706 


14.5364 


13.4113 


lESS 


16.1877 


15.4530 


14.0638 


IRDE 


17.6084 


16.8987 


15.8075 


CIRC 


17.8359 


17.1587 


15.9690 


CMEA 


17.9253 


17.2605 


16.0440 


AEPN 


19.4742 


18.8749 


17.7317 



TabIa 5 Aceleraciones registradas y calculadas del sismo del 12/08/2014 





Distancia 

Epicentral 

(km.) 


Aceleracion 
registrada 

(gal) 


Aceleracion calculada 


Estacion 


C & B 2013 

(gal) 


Zhao et al. 

2006 

(gal) 




PRAM 


13 


48 


60.9 


63.3 


FENY 


12 


38 


64.8 


66.7 


24MA 


15 


18 


54.5 


57.7 


lESS 


16 


31 


51.1 


54.8 


IRDE 


17 


20 


45.9 


50.3 


CIRC 


17 


41 


45.2 


49.6 


CMEA 


17 


31 


44.9 


49.3 


AEPN 


19 


21 


40.4 


45.3 



Antes de comparar las aceleraciones se debe tener en cuenta que la distancia epicentral de cada 
una de las estaciones reportadas por el IPGH no tienen decimales cuando deberian tenerlo debido a que 
se encuentran en diferentes sitios. 

Por otra parte, la base de datos con las que se han obtenido las ecuaciones de C & B y de Zhao 
et al. 2006 presentan para una misma distancia epicentral una gran dispersion de los valores al trabajar 
con el valor medio la probabilidad de excedencia es del 50% si se desea minimizar esta probabilidad se lo 
debe hacer en terminos generales con dos variables denominadas a.s; la primera corresponde a la 
desviacion estandar y la segunda el numero de veces que se considera la desviacion estandar, si £ = 1 o 
2 por ejempio la probabilidad de excedencia se reduce al16 y 5% respectivamente, pero los valores de s 
tambien pueden ser negatives. (Ramirez etal. 2008) Los valores presentados en la tabIa 5 son para e = 0 
con lo que la probabilidad de excedencia es del 50%. 

Las aceleraciones maximas se hallan con el modelo de Zhao et al. 2006 pero son parecidas a las 
que se obtienen con el modelo de C & B 2013, como se aprecia en la figura 8. Se ve que para algunos 
puntos las aceleraciones que reportan los modelos son parecidas a las que se registraron a pesar de que 
en los modelos se considera £ = 0 . 

A 17 km., del epicentre se tienen aceleraciones maximas en campo libre, que varian entre 20 y 41 
gal. Dos de estos registros fueron obtenidos en la parte baja de la ciudad y corresponden a las estaciones 
del Colegio Militar Eloy Alfaro (CMEA) y del Circulo Militar (CIRC) que estan muy cercanas, en estas se 
registro 31 y 41 gal, respectivamente, la restante fue en el IRD ubicado en la Whymper y Coruha (IRDE) 
que se encuentra a mayor elevacion y en mejor suelo. Evidentemente que estas diferencias son por efectos 
de sitio. 
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Figura 8 Aceleraciones maximas encontradas y registradas en el sismo del 12/08/2014. 



Por otra parte, las aceleraciones maximas registradas, a mas de 1 0 km., del epicentre son bastante 
bajas, menores a 0.05 g. Siendo g la aceleracion de la gravedad, que no causan dano en las estructuras 
de hormigon armado y acero que se tienen en la zona donde fueron obtenidos los registros sismicos. 

En la figura 9 se presenta la variacion del PGA (aceleracion maxima) con la distancia Rmp, hallados 
con el modelo de Zhao et al. 2006, para £ = 0, linea contfnua; para £ = 1, linea superior discontinua y la 
£ = -2, linea inferior discontinua, en este ultimo caso la probabilidad de excedencia es del 95%. Se indican 
ademas las aceleraciones registradas. 

1.000 1 1 1 




0.001 



0.000 J 

Rrup (Km) 

Figura 9 Variacion de la aceleracion maxima en funcion de R 
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8. COMENTARIOS Y CONCLUSIONES 



Tres puntos principales han sido presentadas en este articulo y son las aceleraciones maximas 
registradas en campo libre del sismo del 12 de agosto de 2014; los periodos de recurrencia encontrados 
para los segmentos de las fallas ciegas de Quito y el comportamiento de dos modelos de movimientos 
fuertes para determinar: aceleraciones maximas del suelo, formas espectrales y leyes de atenuacion. Con 
relacion a estos puntos se puede indicar lo siguiente. 

Evidentemente en los registros hubo efectos de sitio ya que en lugares que se hallan a la misma 
distancia se tienen aceleraciones muy diferentes en relacion a los valores obtenidos, el valor de la 
aceleracion maxima en la estacion IRDE es menos del 50% del valor registrado en la estacion CIRC, las 
dos estan a la misma distancia epicentral. Queda de manifiesto que los sismos asociados a las fallas 
ciegas de Quito, por tener un foco superficial se atenuan muy rapidamente con la distancia. 

Con relacion a los periodos de recurrencia que son bastante bajos para sismos de gran magnitud 

se puede decir que estos responden al modelo utilizado, en este case fue el del Terremoto Caracteristico, 
que fue ampiiamente descrito en el apartado 4, que no considera la liberacion de energia que se da cuando 
se registran sismos de menor magnitud. 

Finalmente, los modelos de C & B 2013 y el de Zhao et al. 2006 para el presente sismo reportan 
aceleraciones maximas en campo libre, muy similares. Si se incorporaba la incertidumbre en estos 
modelos de seguro que se cubria con las aceleraciones registradas pero el objetivo era comparar los 
resultados medios que reportan las ecuaciones de movimientos fuertes indicadas, con las registradas en 
el sismo del 12 de agosto de 2014. 
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